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I. INTRODUCTION 

The idea of the restoration of spontaneously broken gauge symmetries under extreme conditions (high temperature, 
high fermion densities, strong external gauge fields) was intensely studied in the 1970's and is well known [QJ^] 
(for review see HQ). Under these conditions spontaneously broken gauge theories undergo various phase transitions 
resulting, in most cases, in the restoration of the originally broken symmetries. The consequences of such restoration — 
masslessness of the gauge bosons and fermions, long range character of the interactions, etc. — are very important for 
the properties of the very early Universe. These properties may also play a significant role in the context of neutron 
stars and in the multiparticle production at high-energy interactions. 

It has long been realized that the Lagrangian of the Standard Model with a non-zero vacuum expectation 
value for the scalar field is a covariant generalization of the Ginzburg-Landau equation which describes a wide variety 
of "ordered" systems (ferromagnets, superfluids, superconductors, crystals, etc.) and their phase transitions. From 
this point of view very different systems have common general properties due to their differently revealed broken 
symmetries. Using the analogy with superconductivity, where a strong external electromagnetic current (magnetic 
field) destroys superconductivity resulting in the transition to the normal state, an intense current of weakly interacting 
particles (e.g. of W^-bosons) can restore the gauge symmetry and cause the transition to a new phase with massless 
gauge quanta |||]. 

The present paper is devoted to the study of the restoration of the spontaneously broken gauge symmetry in the 
collisions of intense, energetic W^-boson wave packets. The paper is organized as follows. Section II describes the 
symmetry restoration at zero temperature in the SU(2) Higgs model. In Section III this phenomenon is compared with 
the well-known results of symmetry restoration at high temperature. Section IV is devoted to the stability problem 
of the one-loop radiatively corrected effective potential with the external gauge field in the light of the existence of 
the heavy t-quark. In Section V we present lattice calculations of gauge boson wave packet collisions in which the 
restoration of the SU(2) gauge symmetry is manifested. Section VI contains short concluding remarks. 

II. RESTORATION OF THE BROKEN SU(2) SYMMETRY BY AN INTENSE GAUGE FIELD 

Thispaper is an outgrowth of previous studies of the interactions of gauge field wave packets in the SU(2) Yang- 
Mills H| and Yang-Mills-Higgs Q theories. Being in essense classical, the approach of these publications nonetheless 
reveal interesting phenomena associated with the non-perturbative description of multiparticle production in collisions 
of weakly interacting particles. However, in the present study of symmetry restoration of the scalar field effective 
potential as a result of such wave packet collisions, this treatment is more adequate since the colliding of intense, 
energetic I-U-boson wave packets serve here as a space-time averaged external source exciting the vacuum of the 
Yang-Mills-Higgs system. 

Our consideration is based on the spontaneously broken SU(2) Yang-Mills theory with one Higgs isodoublet coupled 
to the gauge field. This model retains the most relevant ingredients of the electroweak theory. We begin with some 
preliminaries. 

The action describing the model in (3 + 1) dimensions is given by 
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s = 



d 3 xdt 



-\tv{F^F^) + §tr ((£V$)tC$) - A (±tr($t$) 



(1) 



with Dp = 0„ - igA«r a /2, Fpv = F^ v T a /2 = (i/g)[V^V u ], and 

$ = 0° - ir a </> a , (2) 
where r a (a = 1, 2, 3) are Pauli matrices. We work in the unitary gauge where only physical excitations appear: 

$ = (v + p/V2)U(6) , 

(3) 



= U(6)W fi U- 1 (6) - - (8^1/(0)) U' 1 (9) , 

with U{9) — exp(iT a 9 a ) being a local gauge rotation. 

The real field p describes the oscillations of the scalar field about its vacuum expectation value v, and Wp is the 
field of the gauge boson. and p obey the following classical equations of motion: 



[2V FH + M^W V + \fvpW v + \g 2 p 2 W» 



0. 



- /n 1 

(d^ + M%) p + 3V2Xvp 2 + Xp 3 - Z- g 2 vW^W a ^ - -g 2 pW«W a ^ = 



(4) 
(5) 



Dp and F^ v are defined in terms of W M . Mh = 2vvX and M\y = gv/V2 are the masses of Higgs and gauge bosons, 
respectively.]] ^ 

It is easy to see that after a scaling transformation 

x' ll =gvx^ = A^=A, x /v, (6) 

the action (|l|) and the equations of motion (^) and (|J) possess only a single parameter A' = X/g 2 — M\ /8M^. 
However, since in the following simulation of the colliding wave packets there are other parameters involved in the 
initial conditions, we consider the unsealed equations (Q) and (0) in more detail. 

Equation (||) describes the excitations of the scalar field about the Higgs vacuum |$| = v. It has a second exact 
solution p = —\/2v (corresponding to |$| =0) with W® being arbitrary. For this solution the last three terms in (^) 
cancel each other and eq. (Q) describes massless gauge bosons. 

In terms of excitations x = P + v\/2 around |<I>| = 0, the equations (|j) and (||) assume the form: 



M 2 
~~2 



// I , , 9 2 W 2 



8Xv 2 



X + Ax 3 = . 



(7) 
(8) 



The expression g 2 W 2 /8Xv 2 in (|J), where W 2 (x) = W^W a ^ , plays an important role: Depending on the sign and 
magnitude of W 2 (x), the new "ground" state |$| = can be unstable or stable. For positive W 2 (x) this state is 
obviously unstable noting that the effective mass of the scalar excitations described by (|^) is imaginary. 

For time-like (k 2 > 0) gauge bosons W 2 (x) is always negative. In particular, this is so for transverse gauge bosons for 
which Wq = 0. For space-like (k 2 < 0) gauge bosons W 2 (x) can be of both signs. On the other hand, the luminosity 
of transversely polarized gauge bosons in proton-proton collisions at high energy, under realistic conditions, is much 
higher than that of longitudinally polarized ones and increases with energy Q]. Of course, the general suppression by 
a factor M\y/E of the emission amplitude of longitudinal IT-bosons from a quasi-massless fermion is compensated 
by the nonvanishing of this amplitude in the forward direction in contrast to that for transverse M^-bosons leading, 
in general, to the dominance of the cross-sections for longitudinal W-W interactions over transverse ones H. But, as 



1 Note that the gauge field acts as a source for Higgs field excitations in This permits us to consider the VK-field classically 
with respect to the excitation of the scalar field. The equation (^) for gauge field does not possess a source term. 

2 Our normalization of the p-oscillations about v in corresponds to the choice v = 174 GeV. Thus, Gf/V% = l/4w 2 , where 
Gf is the effective four-fermion coupling constant. 
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seen from the above considerations, for the restoration of symmetry the important factor is not the cross section of 
W-W interaction, but rather the intensity of the "pulse" W 2 which appears in equation (||) as an external source (cf. 
the footnote 1). 

Concerning the W-W collisions, in the following we will consider collisions of transverse gauge bosons for which 
the relation d^W? — holds. This also allows us to work in the temporal gauge Wq = 0, which is convenient for 
lattice calculations in the Hamiltonian formulation (see Section V for details). We would like to emphasize that, 
considering the case of transversely polarized W-W scattering, we exclude the important issue of the potentially 
strong longitudinal W-W interaction which may exist beyond the Standard Model ||. The restoration of the broken 
symmetry in the transverse W-W collisions studied here is a phenomenon entirely within the confines of the Standard 
Model. 

Since the colliding wave packets are highly energetic (k 3> M^,M#) and rapidly oscillating in space and time, the 
Higgs field oscillations will only follow effectively the average of the gauge wave packets. We therefore introduce the 
space-time averaged quantity 



Equation (ph describes the excitations of the scalar field x with an effective mass 



V ~ 8\v 2 M 2 - 2 ' 



M 2 H = %V-1). (10) 

For 77 < 1 the excitations around the $ = ground state is tachyonic, indicating instability. For 77 > 1, however, these 
oscillations are stable. In particular, for the spatially homogeneous case there exists an exact stable solution of (g) 

1 SM 2 H X C 

where xo is a constant p0| . 

It is instructive to analyze the effective potential V(x,Tj) for the excitations x described in equation (||): 

V( X ,v) = -^ 2 (l-v)x 2 + ^X i • (12) 

Equation ( |l2| ) has two different stable minima, depending on the magnitude of 77: 

For7?<l, Xmi„ = ±V2v(l - 7y) 1/2 , i.e. |$| = «(1 - ry) 1 / 2 . (13) 
For n>l, Xmi„ = 0, i.e. $ = . (14) 

Stable excitations about these "vacua" have the following squared masses: 

M 2 

M$ = -f\l-v\[l+ 0(1 -f})] , (15) 
Ml r = M* v (l-r 1 )6{l~ V ) . (16) 

Thus for 77 > 1, the broken SU(2) gauge symmetry is restored and oscillations of the scalar field occur about the 
symmetrical state |$| = 0, not about |$| = v. The effective mass of the gauge bosons in the region where the 
symmetry is restored (17 > 1) vanishes. For rj < 1, the vacuum is changed gradually as (1 — rf) 2 , approaching |$| = 0. 
For 77 < 1, the ratio r = Mh/Mw = Mh/Mw has no dependence on 77. 



III. COMPARISON WITH THE HIGH TEMPERATURE CASE 



The above described features are characteristics of a second order phase transition. From this point of view, it is 
illuminating to compare the restoration of symmetry in an intense gauge field associated with time-like VF-bosons 
(k 2 > 0, or W^W afl < 0) with the analogous phenomenon occurring at high temperature. 

For the same SU(2) Higgs model the parameter governing the restoration of symmetry at finite temperature at 
one-loop level is [pl]-^3| 
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v 2 



Ml 
Ml 



(17) 



The high temperature acts here in the same way as an intense time-like gauge field (or as an extra-high fermion 
current density The joint action of these two effects, of course, would make the transition easier. It is important 

to stress that our essentially classical approach is not reliable for rj close to unity where the radiative corrections are 
significant even at small coupling constants. However, it is expected that the radiative corrections are not important 
in the region described by the inequality jl . 



\1-V\ >A, 



(18) 



IV. ON THE ROLE OF THE RADIATIVE CORRECTIONS AT ONE LOOP LEVEL 



Evidence of the existence of the very heavy i-quark (m, t ~ 175 GeV) G3 raises the question of the stability of 
the ground state with and without external influence. In this section we show that the possible inclusion of heavy 
fermions into our radiatively corrected effectived potential (|l^), under realistic conditions outlined in Section 3, does 
not change the resulting picture for moderate rj and the perturbative range of A and g 2 . 

The radiatively corrected potential (|l^) has a well-known contribution from one loop diagrams including the physical 
SU(2) gauge bosons, Higgs scalar, and i-quark |l5|] (we ignore other much lighter quarks and other scalars). We have 
for our choice of the normalization of the x-excitations and of v: 



C x 2 
V c (x, V) = V(X, V) + jX 4 In , 



(19) 



where M is an arbitrary normalization scale which will be chosen in the following equal to v. For C we have U 



C = 



1 



256tt 2 u 4 



4 1 Mjj-AMf 



(20) 



Inserting into (p0|) the physical masses of particles, we are left with 



C = 4 x 1(T 4 



M H 

v 



3.5 



(21) 



As we are interested in the possible destablizing role of the t-quark, we consider the case C < 0, i.e. Mjj < 238 GeV, 
or, A < 0.47. On the other hand, A > 0.03 follows from the present experimental lower bound on the Higgs mass 
(M H > 65.2 GeV) @. 
Now, the minimum of the full effective potential ([11]) is 



v 2 . = 

A mm 



2v 2 (l -rj) 



1 + S 



111 M 2 



0--v) 



(22) 



The term i in the denominator of (^2|), of course, can be absorbed into M. The prefactor of the logarithm in (|2^) is 
small (J^i < 0.05). 

The condition for the validity of the one-loop approximation (we take M ~ v) is: 



Tj > e i c i , 
■KeVIci 



for 1] < 1 
for rj > 1 



(23) 
(24) 



3 The unconventional numerical factor in C is due to our normalization of the p-oscillations in ^ and our choice of v = 174 
GeV. Summation in (|2^) runs over gauge bosons W ± and Z. 
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which is well satisfied considering our restriction ( |18| ) of the whole reliability of the approach and for not very large 
rj. Thus we conclude that the picture does not change drastically by the inclusion of the radiative corrections with 
the dominated /j-quark contribution. 

Indeed, one has for \rj — 1| > A, g 2 and not large rj: 



Mi 



c 



0(1 -n) 



M%=^f\l-V\ 1-t (1 + *(!-»?)) , 



A 

= M^(l - v )[i-£ l n (i _ rj) ) 9(1 - rj) , 



A 



A,= -: 



ld 4 V 



6 dx 4 



= A 



(25) 



From ( |25| ) it is seen that the use of the bare coupling constant A in all previous estimates is reasonable. 

The physical reason for this small influence of the i-quark is clear: the choice of moderate external vector field 
amplitude W 2 (< 0) ensures the reliability of one- loop approximation and, as we see, the radiative corrections cannot 
stop the tendency towards the restoration of the symmetry. 

It is worthy to emphasize that whereas our arguments are well grounded for the rj < 1 region, relations (p5| ) have 
only qualitiative validity even for moderate 77 in the region rj > 1 where infrared singularities connected with the 
masslessness of the gauge bosons are important and practically incurable. These difficulties are common for phase 
transitions in many-body problems treated by perturbation theory or by mean field methods ||. 

Note, finally, that in the potential ( |l9| ) the renormalization scale dependence of its parameters is ignored. However, 
the one-loop potential ( fL9| ) has all the qualitative features of the "full" renormalization group improved potential, and 
the results thus obtained will not be substantially changed if we use the "exact" potential 1 1^] . 



V. RESTORATION OF THE SU(2) GAUGE SYMMETRY IN GAUGE WAVE PACKET COLLISIONS 

As follows from the arguments in Section IV, the reliability of our analysis is confined to the values of not very 
large 77. For rj 3> 1 our treatment is, in principle, not well grounded, but that occurs only in the very early Universe. 
On the other hand, the study of the gradual decrease of the vacuum expectation value as a function of rj from |$| = v 
to |$| = — rj) 1 / 2 for rj < 1 is free from the difficulties mentioned in Section IV. 

In this section we study the collision of gauge wave packets from the point of view of the restoration of the SU(2) 
broken symmetry. In Jjj we have investigated the evolution of the gauge and Higgs's field as a function of the 
colliding wave packet amplitudes and of the mass ratio r = Mh/Ww- For a wide range of the essential parameters 
we have observed collisions of the initial wave packets resulting in final configurations with dramatically different 
momentum distributions. Note that the classical approach used in these calculations, although suffering from the 
common ultraviolet catastrophe of classical field theory revealed in the momentum distribution of the final states (see 
0), is adequate in the present context since we are not interested in the final states. 

Our numerical study here, as in j7|, is based on the Hamiltonian formulation of lattice gauge theory (see [jl9]j|] 
for more details), in which the dynamical variables are link variables defined as 

U e = exp (-igaA c f T c /2) . (26) 

(£ is the link index.) As in J^,|?j we work with a one-dimensional lattice of size L = Na, where N is the number of 
lattice sites and a the lattice spacing. As mentioned in Section II, we "collide" transverse VF-bosons to ensure the 
negativity of W®W afl as well as the temporal gauge condition W§ = 0. 

The initial configuration is given by two well-separated left- and right-moving wave packets originally centered at 
z l(r) with average momenta k = (0,0, k) and width Ak. ^l(_r) are chosen so that the wave packets are positioned 
symmetrically about the center of the lattice: 

^ = (o,o,i,oK#- Zl ,-i) , 

= (0,0,l,0)n c R ii(z - z R ,t) , (27) 
with n£ and being the polarization vectors in isospin space. 
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A right-moving wave packet centered at z — at t = with mean wave number k, mean frequency w = (k 2 + M w ) 1 / 2 
and width Ak is described by 

1 1/2 /"+0O 

t ) = / dk z e -(^-fe) 2 /2(Afc) 2 [cQs( ^ _ kzz)] (28) 



with w = \/k 2 + M w , and 



n = 



(29) 



where the amplitude is fixed by requiring energy equal to fiu per cross-sectional area a. (In the following we set 
h = 1.) From (psj) we obtain 



^| t= o- J^r e - (Afcz)2/2 cos(fcz), (30) 



and 



dr 1 V^(n/w) 



sm(fe) + — — Akzcos(kz) +0 — 



(31) 



Initial condition for the Higgs field is the vacuum solution: 

0° = u, </> a = 0, 0° = 4> a = at i = . (32) 

To determine the number of independent parameters, it is useful to rescale (^) according to (^) in the high-energy 
limit (f2 f=a Co ~ k): 

1,, / 2Afc -{^z'Y/2 fk \ 

^ t/=o = -V» «=o « W ^ r 2 e 9 / cos — z' , (33) 
w V \l^kav l \gv ) 

which contains three dimensionless parameters: k/Mw, Ak/M\y, and aM w /g 0- There appears one more parameter 
in the initial conditions — the relative initial isospin orientation of the colliding wave packets, 9 C . As shown in HQ, 
the parameter # c plays an important role in the collisions: Whereas for 8 C ^ (essentially non-abelian configuration) 
collisions of the wave packets resulted in highly inelastic final states, for parallel isospins (pure abelian configuration) 
the collisions were practically elastic. However, for the phenomenon we are interested in here — restoration of the 
broken symmetry — the coupling between the gauge and Higgs field is essential and not the relative orientation of the 
gauge bosons in isospin space. 

It is convenient to have a scale parameter in the dimensionless parameters v rather than Mw- Combining equations 
of motion and initial conditions, we finally have four independent parameters: 

r.£L,i,^ lM d2£. (34, 
Mw v v g z 

In our simulations we always chose k 3> Mw and k 3> Ak to model high energy collision and to have wave packets as 
well-defined objects. The gauge coupling constant was fixed to g — 0.65. 

We estimated the parameter r\ using two different configurations of the gauge boson wave packets. Firstly, we 
consider the co nfig uration of the two incident wave packets at the space-time point when they overlap completely. 
Then from (||),(^), and we have 

p-K 2 , for parallel isospin orientation 

^k 2 , for orthogonal isospin orientation (35) 
, for antiparallcl orientation 

where just the parameter we have introduced earlier (see ( |32"| ) in ||). As found there, for k > 1 the 

W-W wave packet collisions produced strongly inelastic final states (for non-parallel isospin orientation). However, 
one should note that according to Eq. (||) the scalar field reacts to the gauge field long before the overlapping of 



G 



the wave packets (see also figures below). That makes the behavior of the collisions less sensitive to the relative 
orientations of the wave packets in isospin space. Practically, as seen in our simulations, there is no difference between 
the parallel and the antiparallel cases. 

Let us estimate 77 taking the initial (t — 0) configurations of two separate wave packets. From (§),(||), and (|J) we 
have after spatial averaging 

K ' 2 2Afc 

V=~T= r— .9,9 ■ ( 36 ) 



Note that no isospin factor appears here after spatial averaging. As we see, the two estimates (|35|) and (|36|) almost 
coincide. A sufficiently large 77 is achievable for realistic values of r > 1 and permissible ratio Ak/k if we take a < 1/v 2 . 

The natural spread of the W^-boson wave packets emitted by fast moving fermions (quarks, leptons) is of order 
v in the comoving reference frame. Therefore, the transverse area of the VF-boson wave packet is not larger than 
1/v 2 . Of course, the more reliable estimate of 77 depends on the detailed shape of the wave packet and requires a full 
three-dimensional analysis. 

Figures 1 and 2 describe the results of our lattice simulations of the collisions of W wave packets. In view of the 
above mentioned relation between the parameters 77 and K (or k') (see ( |35|) and (^6|)) we present here only results 
for wave packets with orthogonal isospins since parallel isospin orientation leads qualitatively to the same symmetry 
restoration patterns. 

In the first column of Figure 1, the space development of the colliding W wave packets with 77 = 1.32 is shown for 
different times. As seen from this figure, at time t « 300 the W wave packets collide and then begin to separate. 
Just about at this time one expects to observe the restoration of symmetry, i.e. the oscillations of the scalar field 
|$| about a new ground state located below the "old" vacuum \&\/v = 1. After the separation of the wave packets 
(t > 300) the excitations of the scalar field tend again to oscillate about the "old" vacuum, i.e. the gauge symmetry 
is broken again. This is seen clearly in the second and third columns of Figure 1. The second column shows the 
space-time evolution of the magnitude of the Higgs field |3>|/w. The third column shows the Higgs field smoothed 
over 50 lattice sites, suppressing the small scale fluctuations, in order to facilitate a comparison with the definition 
@ of the effective parameter 77 in terms of the averaged strength of the VK-boson field. At the moment of interaction, 
t « 300, the lowering of the average value of the Higgs field at the location of the coalescing wave packets is clearly 
visible. 

In Figure 2, we show the 77 dependence of the W-W collisions (first column) and the corresponding behavior of the 
scalar field (second and third columns) at t = 300. To change 77 we vary the amplitude of wave packets by changing 
a; all other parameters are kept fixed. It is evident that no significant symmetry restoration occurs for 77 < 1, as 
expected. For 77 3> 1 the situation becomes more complicated, probably because already a single wave packet can 
restore the symmetry in this case. It is seen, however, that the spatial region with restored symmetry becomes wider 
as the parameter 77 grows. 



VI. CONCLUDING REMARKS 



The above study of colliding gauge wave packets in the spontaneously broken SU(2) Yang-Mills-Higgs model visu- 
alizes the process of the SU(2) gauge symmetry restoration in some finite space-time region as a result of the presence 
of the intense gauge pulses. In some sense, it qualitatively resembles the situation realized in the very early Universe 
with the Big Bang scenario where the restoration is caused by the high temperature. The approach used in this paper 
might be applied to the study of chiral symmetry restoration. As is well-known, there exists an isomorphism between 
the Higgs sector of the standard model and the linear cr-model. This isomorphism is realized through the following 
identifications: Wl <-> 7T, H <-> a, v «-> /„•, where Wl is the longitudinal gauge field and H the physical Higgs field. 
fn is the 7r-decay constant defining the scale of the chiral dynamics. As a consequence of this isomorphism, one might 
expect that the spontaneously broken chiral symmetry can also be restored in the presence of a strong external pion 
field and such a physical setting might be found in heavy-ion collisions or in the interior of nuclei and neutron stars. 

We have seen in Section V that the restoration of the gauge symmetry in W^-wave packet collisions occurs at 
amplitudes for which strongly inelastic behavior has also been observed. To achieve a sufficient value of the amplitude 
is not easy in practice. The problem, presumably, is connected with the general problem of the classical (quasi- 
classical) treatment of the amplitude for the transition from a few-particle to a many-particle final state (see, e.g. 
pfj|). If, as conjectured in |H| (see also p2]-|2^]), there exist classical trajectories interpolating between incoming 
few-particle and outgoing multi-particle states, our results would lead us to expect that these trajectories will also 
exhibit the phenomenon of symmetry restoration in the transition region. It is also worthwhile mentioning that the 
classical approach is more adequate in the context discussed here, where the colliding intense, energetic VF-boson 
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wave packets serve as a classical external source exciting the ground state of the Yang-Mills-Higgs system, and we are 
not interested in the specific nature of the final states. 

In conclusion, we emphasize that the restoration of the gauge symmetry discussed above does not depend on 
the dimensionality of space. It is expected that the symmetry restoration would not persist as long in three space 
dimensions because the wave packets disperse more rapidly after the collisions, causing the average squared amplitude 
((W"") 2 } to decay more rapidly. The same arguments can apply to the problem of multiparticlc production in wave 
packet collisions [|]|7]]. This issue, as well as the problem of the reliability of the classical treatment of the (few— > 
many) particle amplitude makes the study of the gauge wave packet collisions in (3 + 1) dimensions desirable. 
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FIG. 1. Space-time development of symmetry restoration induced by two colliding gauge wave packets with orthogonal 
isospins in the presence of the Higgs vacuum condensate. First column shows the scaled gauge field A\/v as a function of space 
coordinate z at five chosen times. Second column demonstrates the corresponding space-time evolution of the scaled Higgs 
field |$|/f. Third column shows the scaled Higgs field after smoothing over 50 lattice sites. This simulation as well as the 
following (Fig. 2) was done on a lattice of sites iV = 2048 and lattice spacing a = 1. The parameters were k — n/4, Ak — tv/16, 
M h = M w = 0.15, g = 0.65, and a = 1. 
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FIG. 2. j)- dependence of the symmetry restoration. We demonstrate at a chosen time t = 300 the space pictures of the 
colliding gauge wave packets and the Higgs field for four different values of 77. First, second, and third columns shows the scaled 
gauge field |A|/v, the scaled Higgs field and the smoothed Higgs field respectively, each as a function of the space 

coordinate 2. All parameters were the same as in Fig. 1 except a, which was to varied to arrive at different values of r\. 
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